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^ Abstract 

The optimal diversity-multiplexing-delay tradeoff for the multi-input multi-output (MIMO) automatic repeat request 
(ARQ) channel can be achieved using incremental redundancy lattice space-time codes coupled with a list decoder for 
joint error detection and correction. Such a decoder is based on the minimum mean-square error lattice decoding principle 
which is implemented using sphere decoding algorithms. However, sphere decoders suffer from high computational 
'j- ^ ' complexity for low-to-moderate signal-to-noise ratios, especially for large signal dimensions. In this paper, we would 
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'"1 Uke to construct a more efficient decoder that is capable of achieving the optimal tradeoff with much lower complexity. In 
^ O ^ particular, we will study the throughput-performance-complexity tradeoffs in sequential decoding algorithms and the effect 
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of preprocessing and termination strategies. We show, analytically and via simulation, that using the lattice sequential 
decoder that implements a time-out algorithm for joint error detection and correction, the optimal tradeoff of the MIMO 
ARQ channel can be achieved with significant reduction in decoding complexity. 

I. Introduction 



^ Automatic Repeat reQuest (ARQ) is an efficient communication strategy that uses feedback to achieve high reliability, 

j> and is widely used in many wireless networks (e.g., LTE and WiMAX) (refer to [IJ for a detailed study about several 
• ^ 

^ ARQ schemes). In its early stages, ARQ was used in conjunction with codes with good error detection capabilities. 
However, such codes increase the number of retransmissions which significantly reduce transmission rate (throughput) and 
increase delay. This may become undesirable for many communication systems, particularly in wireless fading channels. 
To overcome such problems, hybrid-ARQ system was introduced which uses forward error correction techniques (e.g., 
block and convolutional codes) |l2]|-|l7]|. This, however, comes at the expense of increasing the complexity of the receiver 
The design of low complexity receivers for ARQ systems that achieve near optimal performance and high throughput 
is considered a challenging problem. 



The class of sequential decoders is among the most promising decoders that can handle high data rates with low 
decoding complexity. Sequential decoders lEl, ||9l are tree search algorithms that were originally constructed to decode 
convolutional codes transmitted via discrete memoryless channel, and are well-known to achieve near optimal performance 
with very low (average) decoding complexity. However, there is still a non-zero probability that the decoding complexity 
(time) becomes excessive, especially when the channel is very noisy. In this situation, the decoder encounters buffer 
overflow that results in a decoding failure. It is this probability that limits the performance of the sequential decoder. 
Fortunately, the decoding failure probability can be totally eliminated using systems with feedback channel. For that 
reason, sequential decoders were adopted with ARQ systems iBj-lHl due to their ability to detect for retransmission before 
ending the decoding search which results in huge saving in decoding complexity while maintaining high throughput. All 
of this makes sequential decoding very promising and attractive for use in systems with repeat request. 

Many sequential decoding algorithms (e.g., stack algorithm ||9l) were modified for the use of signal detection and 
decoding in ARQ systems. Among those algorithms that is considered simple but efficient is the so-called time-out 
sequential decoding. In this algorithm, the decoder simply tracks the number of computations performed by the decoder 
and asks for retransmission if the computations become excessive and exceed a certain predetermined time limit. This 
results in reducing the decoding complexity by terminating the search during high channel noise. For the case of single- 
input single-output discrete memoryless channel, it was shown (see lH) that there exists an optimal time-out limit value 
that maximizes both performance and throughput while achieving low decoding complexity. In this paper, we would 
like to extend the work in [3| to the multi-input multi-output (MIMO) ARQ channel, particularly, to the quasi-static, 
Rayleigh-fading MIMO ARQ channel. In particular, we will study the throughput-performance-complexity tradeoffs in 
sequential decoding algorithms and the effect of preprocessing and termination strategies such as the time-out algorithm. 

Diversity-multiplexing trade-off (a result that has been developed by Zheng and Tse [10] which shows a rigorous 
fundamental tradeoff between the data rate increase possible via multiplexing versus the channel error probability reduction 
possible via diversity) has become a standard tool to evaluate the asymptotic performance of the MIMO channels at 
the high signal-to-noise ratio (SNR) regime. Space-time coding, a powerful coding technique, is used to exploit such 
tradeoff ifTTI -lfT?!. In MIMO ARQ systems, the delay introduced by the channel provides a third dimension in the 
tradeoff region. To achieve the optimal diversity -multiplexing-delay tradeoff of the MIMO ARQ channel. El Gamal et. 
al. |[T5l proposed an efficient coding scheme called Incremental Redundancy LAttice Space-Time (IR-LAST) coding. 
These lattice-based construction of space-time codes were designed using linear random coding techniques | fT4l . The 
problem of constructing explicit optimal IR-LAST codes for the above mentioned MIMO ARQ channel was discussed 
in |[T6l . However, in both papers, a list lattice decoder for joint error detection and correction (usually implemented via 
sphere decoding algorithms lITTl ) is an essential part for achieving the optimal tradeoff. The draw back of using such 
a decoder is that at low-to-moderate SNR and for large signal dimensions, the size of the candidate list could become 



extensively large. This motivates us to search for a more efficient joint decoding technique that is capable of achieving 
the optimal tradeoff with a fairly low decoding complexity. 

In this paper, we propose a different approach for joint detection and decoding based on the efficient lattice stack 
sequential decoder f9l. Such a decoder was developed in f^Ol for solving the closest lattice point search problem, which 
is related to the optimum decoding rule in MIMO channels. We implement the time-out algorithm at the decoder to 
predict in advance the occurrence of high channel noise. This results in less wasted time trying to decode a noisy signal 
and hence improving upon decoding complexity. We show, analytically and via simulation, that the optimal tradeoff can 
be achieved using such a decoder with significant reduction in (average) decoding complexity. 

The organization of this paper is as follows. In Section II we introduce our system model and review basic concepts of 
lattices and the IR-LAST coding scheme. In Section III, we briefly describe the mechanism of lattice sequential decoder 
and derive lower and upper bounds for the performance of such decoder from a lattice point of view. These bounds are 
considered the primary elements in constructing our new MIMO ARQ decoder. In section IV, we introduce the time-out 
algorithm and prove its optimality in terms of the diversity-multiplexing-delay tradeoff. Finally, we verify our theoretical 
results using simulation in section V. 

Throughout the paper, we use the following notation. The superscript denotes complex quantities, ^ denotes transpose, 
and ^ denotes Hermitian transpose. We refer to g{z) = as lim^_i.oo 9{z) / log{z) = a, > and < are used similarly. For 
a bounded Jordan-measurable region TZ c M"^, V{TZ) denotes the volume of TZ. We denote 5m(r) by the m-dimensional 
hypersphere of radius r with V{Sm{r)) = (7rr^)™/^/r(m/2 + 1), where T{x) denotes the Gamma function. Also, 1^ 
denotes the m x m identity matrix, and (g) denotes the Kronecker product. The complement of a set A is denoted by A. 



A. ARQ MIMO Channel 

Consider a MIMO ARQ system with M-transmit and A'^-receive antennas, a maximum of L rounds, no channel state 
information (CSI) at the transmitter and perfect CSI at the receiver. For the MIMO ARQ channel model, we follow in the 
footsteps of El Gamal et al. fTSl and use the incremental-redundancy ARQ transmission scheme. We restrict ourselves 
to the one-bit feedback (ACK/NACK) MIMO ARQ model. The ARQ feedback channel is assumed to be zero-delay and 
error-free. The complex baseband model of the received signal at the ^-th round can be mathematically described as 



(SNR) per receive antenna. The elements of both the noise matrix and the channel fading gain matrix are assumed to 



II. System Model 




(1) 



where G ^MxT jj^g transmitted signal matrix, T is the number of channel uses, Y'j, G C^^^ is the received signal 
matrix, VF^ G ^NxT jj^g noise matrix, H^^ G £^NxM jj^g channel matrix, and p is the average signal-to-noise ratio 



be independent identically distributed zero-mean circularly symmetric complex Gaussian random variables with variance 

In this paper, we assume two different scenarios of channel dynamics. The first model being the long-term static 
channel, where the channel coefficients remain constant during all L rounds, i.e., ff^ = H'^ for all 1 < ^ < L. The 
second scenario is the short-term static channel, where the channel remains constant during each round and changes 
independently at each round. Also, the following short-term average power constraint on the transmitted signal is assumed 

IE{ IIX^Ill.} < MT. (2) 

The equivalent real-valued channel model, after I transmission rounds, corresponding to ([T]) can be written as 

yi = Hix + Wi, (3) 

where we define x = (x^.^, . . . ,xj ^ . . . ,a;J j,)""", with xj^ = {R{[Xj]t}'^ ,Q{[Xj]t}'^y , and w = {wj^^, . . . ,wj^^, 
. . . ,wjrp)'^ , with wj^ = (R{[Wf\t}'^ , Q{[Wf\t}'^)'^ . The vector G M^^^^ represents the total signal received over all 
transmitted blocks from 1 to i. The equivalent real-valued channel matrix, H£, has dimension 2NT£ x 2MTL and is 
formed by taking the first 2NTi rows of the matrix which is composed by L diagonal blocks, where each diagonal 
block takes the form 

There are several interesting issues for why we resort to the real- valued channel model. First, the input-output relation 
describing the channel that is given in Q allows for the use of lattice theory ll26l which may simplify our analysis. Second, 
space-time codes based on lattices have been widely used in MIMO channels due to their low encoding complexity (e.g., 
nested or Voronoi codes) and the capability of achieving excellent error performance |14|. Another important aspect of 
lattice space-time (LAST) codes is that they can be decoded by a class of efficient decoders known as lattice decoders 
(implemented using sphere decoding algorithms). 

B. IR-LAST Coding Scheme 

An m-dimensional lattice code C{A,Uo,Tl) is the finite subset of the lattice translate A + «o inside the shaping region 
TZ, i.e., C = {A + Uo} n7^, where 7^ is a bounded measurable region of M™. It is well-known lO that an (M xT)x L 
space-time coding scheme is a full-dimensional LAST code if its vectorized (real) codebook (corresponding to the 
channel model ^) is a lattice code with dimension m = 2MTL. Let Ac = {x = Gz : z € Z™} be a lattice in 
where G is an m x m full-rank lattice generator matrix. Denote (5a(x) = argmin_xgyv 1-^ ~ x\ the nearest neighbour 
quantizer associated with a lattice A. The Voronoi cell that corresponds to the lattice point a; € Ac is the set of points in 




closest to X, i.e., V^(G) = {r € M"" : QA^r) = x} and its volume is given by Vc = V{Vx{G)) = y/detiC^G). The 
covering radius rcov(G) is the radius of the smallest sphere centered at the origin that contains Vo(G). The packing radius 
^'packCG) of the lattice Ac is the radius of the largest sphere centred at the origin that is within Vo{G). The effective radius 
refr(G) is the radius of the sphere with volume equal to Vo{G). Let be a sublattice of Ac with fundamental Voronoi 
region (corresponds to lattice point 0) V^. An m-dimensional nested lattice code is defined by C = {Ac + Uq} nVs- 

We say that a LAST code is nested if the underlying lattice code is nested. Here, the information message is effectively 
encoded into the cosets A^ in Ac. As defined in [14], we shall call such codes the mod- A scheme. The proposed mod- A 
scheme works as follows. Consider the nested LAST code C defined by Ac (the coding lattice) and by its sublattice Ag 
(the shaping lattice) in M*". Assume that A^ has a second-order moment (7^(As) = 1/2 (so that it uniformly distributed 
over Vs satisfies E{|Mp} = MTL). The transmitter selects a codeword c G C, generates a dither signaQu with uniform 
distribution over V^, and computes x = [c — u] mod A^. 

For the MIMO ARQ channel, we use the mod- A incremental redundancy scheme that was provided in |[T5l . The 
signal X is partitioned into L vectors of size 2MT each. Those vectors are transmitted, sequentially, in the different 
ARQ rounds based on the ACK/NACK feedback. Upon completion of the £ < L transmission, the receiver attempts to 
decode the message using lattice stack sequential decoder (see Section III) that implements a sort of deadline algorithm. 
In particular, the received signal, y^, is multiplied by the forward filter matrix Fi of the minimum mean-square error 
decision feedback equalization (MMSE-DFE) corresponding to the truncated matrix Hg, and then the dither signal filtered 
by the upper triangular feedback filter matrix Bi of the MMSE-DFE is added to it (the definitions and some useful 
properties of the MMSE-DFE matrices F, B are given in |[T4ll ). In this case, the received signal can be expressed as 

y'i = Feye + Beu = Bid + e', (4) 

where c' = c + A, A = — Qa, (c — «), and e' = —[Bi — FiHi]x + FiWi. One can easily verify (see |[l4l) the relationship 
between and the channel matrix through the following equations: 
For the case of long-term static channel, we have 

det(BjB,) = (det (/ + j-^{H')^H-^y , (5) 

and for the short-term static channel, we have 

£ 

detiBjBe) = H det (/ + -^{H^)'' H'^^^'^ . (6) 

The basic idea in this approach is to use a modified lattice stack sequential decoder for joint error detection and 

'a dither signal is a random signal that is used to make the MMSE estimation error independent of the transmitted codeword (see |18| for 
further details). 



correction. The decoder first check if the channel is in outage. In this case, an error is declared and a NACK is sent 
back. If not, we use the modified lattice sequential decoder to find a lattice point that satisfies a certain predetermined 
condition (e.g., a time-out limit). Now if no point is found, an error is declared, and hence, a NACK bit is fed back. If 
a point is found to satisfy such condition then we proceed to the next step to find the codeword as c = [Gz] mod A^. 
The only exception to this rule is at the L-th ARQ round, where the regular lattice stack sequential decoder is used to 
find the closest lattice point. 

C. Rate and Reliability 

Let 7] be defined as the average throughput of the ARQ scheme, expressed in transmitted bits per channel use. Following 
the definition used in [15], rj can be expressed as 

where Ri denotes the rate of the first block in bits per channel use, p{€) = Pr(^i, . . . , Ai) with Ag, denoting the event 
that an ACK is fed back at round I. 

Let El denotes the event that the transmitted message is incorrectly decoded by the ARQ decoder, then the probability 
of error can be upper bounded as 

L-l 

Pe<Y. + Pr(^L), (8) 

1=1 

where Pr(ii^^,^^) takes the definition of the probability of undetected error at round I < L — \. 
Define respectively the effective ARQ multiplexing gain and effective ARQ diversity gain as 

p-5>oo log p p-s>oo log p 

Achieving the optimal diversity-multiplexing-delay tradeoff in ifTSl was performed using an incremental redundancy 
ARQ scheme coupled with a list lattice decoder for joint error detection and correction. This decoder (corresponds to Q) 
finds all lattice points that satisfy (see [15]) 

{a; gM™ : \y' - Bexf < MTL{1 + -flog{p))} , 

where 7 is a constant chosen appropriately to ensure the achievability of the optimal tradeoff. This decoder can 
be efficiently implemented using sphere decoding algorithms (see for example iflTl ). Sphere decoders, however, are 
computationally very complex especially for low-to-moderate SNR and large signal dimensions where the output of the 
list sphere decoder can become extensively large that may result in a waste of time trying to decode the message. Hence, 
it is of great interest to search for a low complexity joint detector and decoder that can achieve the optimal tradeoff. This 
fact motivates us to replace the list lattice decoder by a more efficient retransmission strategy using lattice sequential 



decoders Il20l . Il2ll . The strategy is based on the stack algorithm and is designed to predict the occurrence of an error 
in advance by monitoring the number of computations performed by the decoder. Before we do that, we would like to 
introduce next the sequential decoder for lattice codes and some of its parameters that will be used to construct our new 
MIMO ARQ joint error detection and decoding scheme. 

III. Lattice Sequential Decoder: Performance Bounds and Complexity Distribution 

A. Lattice Stack Algorithm 

The Stack sequential decoder is an efficient tree search algorithm that attempts to find a "best fit" with the received 
noisy signal. As in conventional stack decoder ||9l, to determine a best fit (path), values are assigned to each node on the 
tree. This value is called the metric. For the lattice stack sequential decoder, the metric at the ^-th round [corresponds 
to Q] is given by (see |[20l ) 

^,{zle) = bk-\y';'l-Rl2z'^\^ yi< k<m, (9) 

where = [z^, - ■ ■ ,Z2,zi]'^ denotes the last k components of the integer vector z € Z™, il^^j is the lower k x k matrix 
of Ri that corresponds to the QR decomposition of the code-channel matrix B^G = QeRe at the i-th round, y"^ is the 
last k components of the vector y^' = Qjy'^, and 6 > is the bias term. The bias parameter is critical for controlUng 
the amount of computations required at the decoding stage. 

As the decoder searches nodes, an ordered list of previously examined paths of different lengths is kept in storage. 
Each stack entry contains a path along with its metric. Each decoding step consists of extending the top (best) path in 
the stack. The determination of the best and next best nodes is simplified in the closest lattice point search problem by 
using the Schnnor-Euchner enumeration 11171 which generates nodes with metrics in ascending order given any node z'f. 
The decoding algorithm terminates when the top path in the stack reaches the end of the tree (refer to ||9l and |[20l for 
more details about the algorithm). 

The main role of the bias term b used in the algorithm is to control the amount of computations performed by 
the decoder. In this work, we define the computational complexity of the joint lattice sequential decoder as the total 
number of nodes visited by the decoder during the search, accumulated over all ARQ rounds, until a new transmission 
is started. Also, the bias term is responsible for the excellent performance-throughput-complexity tradeoff achieved by 
such decoding scheme. The role that the bias parameter plays in the new efficient decoding algorithm will be discussed 
in details in the subsequent sections. 

B. Performance Analysis: Lower and Upper Bounds 

In this section, we would like to derive lower and upper bounds on the error performance of the lattice sequential 
decoder. These bounds work as the primary elements for constructing our new efficient decoder for the MIMO ARQ 
channel. 



Consider the detection at the ^-th ARQ round. For simplicity, we consider here the long-term static channel (similar 
arguments can be done for the short-term static channel). Assume the received signal is = B^x + ee, and denote 
EidiBg) and Esd{Bi,b) as the events that lattice decoder and lattice sequential decoder make an erroneous detection, 
respectively, where b is the bias term that was introduced in (|9]). It is well-known |[2ll that lattice decoders outperform 
lattice sequential decoders for any b > 0. In fact, the performance of the lattice decoder serves as a lower bound of the 
lattice stack sequential decoder. Now, lattice decoding disregard the boundaries of the lattice code and find the point of 
the underlying (infinite) lattice closest to the received point. As such, due to lattice symmetry, one can assume that the 
all-zero lattice point is transmitted. For a given lattice Ac, we have 



P{Eia{Be)\Ac) = Pr ( J {2{B,xyei > \Bex\^} j < P{EUBe,b)\A,), 



(10) 



where A* = Ac\{0}. 

For the lattice sequential decoder, it seems a bit difficult to obtain an exact expression for its decoding error probability. 
Instead, we seek to derive an upper bound for the error performance of such a decoder which can be done as follows: 

P{E,d{Be, 5) I Ac) < Pr j \J {/i(z, £) > /i,ninW} 

< Pr j U {\Bix\^ - 2iBex)'^ee < bm} | (11) 



Pr J hiBex^ee > \Bex\^ ( 1 




where (a) is due to the fact that in general, fi{z,l) > ^min(^) is just a necessary condition for x = Gz to be decoded 
by the stack decoder, ^min = min{0, b — |e^Jp, 26 — je^^P, • • • , bm — |e^™p} is the minimum metric that corresponds 
to the transmitted path with e'^ = Q^ee, and (6) follows by noticing that — (/Xmin + I^^P) < 0. Following the footsteps 
of im, one can show that ([TT]) can be finally upper bounded as 



P{Esd{B^, b)\Kc) < P{EidiBi)\A,), (12) 

where 

Bi, (13) 



where -Rmod(^) is the rate at round i that can be achieved using MMSE-DFE lattice decoding and according to ^ is 



given by 



R is the transmission rate, and = 0.5(2rpack(-BfG)/rcff (.B^G))^ . Interestingly, one may show that is lower 
bounded by a constant independent of SNR and I and as a result it has no effect on the performance in the SNR scale 
of interest. 

It is clear from the above analysis that the lattice stack sequential decoder approaches the performance of the lattice 



decoder as 6 — )• 0, i.e., Esii{Bi,0) = Eia{Bg). Moreover, the upper bound (12) corresponds to the probability of 



decoding error of a received signal ye = B^x + decoded using lattice decoding and is valid for all values of 
5 < 2[^»°<i(^)--R]/A/0(-£-)^ j E^^(Bi,b) = Eid{Bi). Therefore, for a given lattice Ac, channel matrix H^, and a bias 
term 6 > 0, one can bound the error performance of the lattice sequential decoder as 



P{Eid{Be)\A,) < P{Esd{BtM^c) < P{Eid{Be)\K). 



(14) 



By averaging (14i over the ensemble of random lattices Ac, one can show that (see [21 j) for a fixed non-random 
channel matrix ff^, the rate 



Rb{Hl,p) ^ max|ii^od(^^^p) - 2MLlog (l + ^l+M^^^ ,o| 



(15) 



is achievable by LAST coding and MMSE-DFE lattice Fano/Stack sequential decoding with bias term & > 0, where a 
is given by 



a(£) 



rosjBeG) 
2rpack(5^G) 



h. 



(16) 



The equations ( 15 1 and ( 16 1 suggest that as long as the channel is well-conditioned, one may use large of values of bias 
term which is needed to achieve low decoding complexity (as will be discussed later). On the other hand, if the channel 
is close to outage, very low values of b must be chosen in order to maintain high achievable rates. For example, if 6 = 0, 



(15 1 reduces simply to Ri, = R^od^ which corresponds to the rate achievable by the MMSE-DFE lattice decoder. 

In general, for the long-term static channel, following the footsteps of [21.1 . one can show that if h is allowed to vary 
with SNR and the channel statistics as 

/ X t/2ML 

1 - 



infii(i + A)^/^^ 



M 

n(i+pAo 

i=l 



f 2rp,,^{BeGy ^ 
V r,s{BiG) 



(17) 



In this case, one can easily show that by substituting 6 in ( 15 1 and ( 16 1, we get 



Rb{x,p) = nog7]ix,p). 



(18) 



The term 7]{\, p) can be chosen freely between 1 and Ili^iCl+pAi) ^^^^ maximum achievable rate under lattice decoding). 



Depending on the value of ri{X, p) we obtain different achievable rates and hence different outage performances. 

We define the outage event under lattice sequential decoding as Ob{p, i) = {Hf : Rb{H'^, p) < R}. Denote R = r log p. 
The probability that the channel is in outage, Pout(^, b) = Pr(C'b(/>, £)), can be evaluated as follows: 

Pont{i,b) = Pr(^logr7(A,p) < R) = p-Mr/i)_ (19) 

where di,{r) is the outage SNR exponent that is achieved by the MIMO channel with no ARQ (i.e., with L = 1). For 
simplicity, we may express 

M 

r/(A,p) = </.n(l + pA,)«% (20) 
1=1 

where < < 1 is a constant independent of p, and Ci) ^1 ^ ^ ^ are constants that satisfy the following two 
constraints: J2i=i d — ^^'^ Ci ^ C2 > • • • > Cm > 0. It is not a simple matter to find a closed expression of dh{r) 
for all possible values of Q. However, one case of special interest, where a closed form of db{r) can be easily derived, 
corresponds to d G (0, M). In this case, one can show [see EH, Theorem 4] that for the MIMO channel with no ARQ, 
db{r) is the piecewise-linear function connecting the points {r{k),d{k)), k = 0,1, - ■ ■ , M, where 

M 

r(0) = 0, r{k) = ^ Ci, 1 < ^ < 

i=M-fe+l (21) 

d{k) = {M - k){N - k), < k < M. 
The above choice of Q leads to a bias term (assuming high SNR) 

&(A,p)-^n(l + pA.)(^-^-)^/*^^ (22) 

i=l 

which may grow exponentially with SNR as p'^, e > 0, when the channel is not in outage. This may significantly 
reduce the decoding complexity. However, although the diversity at r = is not affected by the coefficients Q 
{db{0) = MN), the more unbalanced the coefficients are, the worse the DMT is. It is important to note that, other than 
the uniform assignment Ci = C2 = • • • = Cm = 1> the optimal tradeoff 

d*{r) = {M -r){N -r), V r G [0, min{M, iV}), (23) 

cannot be achieved. In this case, we have h = — 0^/^*^], i.e., 6 is a constanj^ independent of p. In other 

words, one cannot let b to scale with SNR as p'^ (to achieve very low decoding complexity) if the optimal DMT is to 
be achieved. Therefore, by varying the decoder parameter (bias term), one gets different performance-rate-complexity 
tradeoffs. 

^For the short-term static channel, although the channel changes from round to round, achieving the optimal tradeoff requires also the use of 
a constant bias term that is independent of the round I. Otherwise, for non-fixed b, the bias term must be adapted with the channel state at each 
ARQ round. 



It is a simple matter to extend the above result to the i-th ARQ round. One can show that, there exists a sequence of 
full-dimensional LAST codes with block length T > (M + N — l)/£ that achieves the DMT curve dh{ri/i) which is 
the piecewise-linear function connecting the points {r{k),d{k/i)), k = 0,1, - ■ ■ , M, where 

M 

r(0) = 0, r{k) = ^ Ci, 1 < ^ < Af, 

i=M-k+i (24) 
k\ ( k^ 



d{klt) =\^-j^j\N-jy 0<k<M. 
In this case, one can show that there exists a lattice code Ac such that 

P{Esd{Be,b)) = p-''^('^/'\ yi<£<L. (25) 

under the condition > + M - 1. 

It must be noted that, except for the case of i = L (the final ARQ round), db{ri/i) is not the best achievable SNR 
exponent for the MIMO ARQ channel at any round £ < L. This is because in the above analysis we have ignored the 
detection capabilities that the MIMO ARQ can have. By carefully designing the lattice sequential decoder to include 
an efficient error detection mechanism, we will show how the SNR exponent can be raised up to achieve the optimal 
tradeoff of the channel at any round £. 

From a lattice point of view, the event of error under lattice decoding can be expressed as the event that the received 
signal is located outside the fundamental Voronoi region of the underlying (infinite) lattice. In this case, assuming was 



transmitted, one may express the bounds in ( 14 1 



P{e i Vo(S,G)|A,) < P(^,rf(B,,6)|Ac) < P(e i Vo(S,G)|Ac). (26) 



where Bi is as defined in ( 13 1. Interestingly, the upper bound in (26 1 provides us with the fact that as long as e G Vq{BiG), 
the received signal is correctly decoded using lattice sequential decoding. As will be shown in the sequel, this fact can 
be used as the basic tool to construct our new efficient joint detection and decoding scheme. 

C. Computational Complexity Distribution 

An important parameter of the lattice sequential decoder is the distribution of computation, which characterizes the time 
needed to decode a message. It is well-known ||9l that the number of computations required to decode a message using 
sequential decoders is highly variable and assume very large values during intervals of high channel noise. Moreover, 
due to the random nature of the channel matrix and the additive noise, the computational complexity of such decoder 
is considered difficult to analyze in general. However, in ll2ll the authors have shown that, for the MIMO channel with 
no ARQ under lattice sequential decoding with constant bias term, the tail distribution becomes upper bounded by the 
asymptotic outage probability with SNR exponent that is equivalent to the optimal diversity-multiplexing tradeoff of 



the channel. This upper bound is shown to be achieved only of the number of computations performed by the decoder 
exceeds a certain limit when the channel is not in outage. This result can be easily extended to the MIMO ARQ channel 
as will be explained in the following. 

Consider again decoding the received signal at the ARQ round I = 1, • • ■ , L. Let (j){z\^t) be the indicator function 
defined by 

1, if node z\ is extended; 



(27) 
0, otherwise, 



and let Mj{t) be a random variable that denotes the total number of visited nodes during the search up to dimension j, 
at round ^ < L. In this case, Nj {£) can be expressed as 

^^■W = E E '^(^^^)- (28) 

k=l z'leZ'' 

For the case of MIMO channel with no ARQ (i.e., L = 1), following the footsteps of [21], one can show that 
the asymptotic computational complexity distribution of the MMSE-DFE lattice sequential decoder (assuming no error 
detection) is given by (for T > N + M - 1) 

> r) < Pr(AA^ > r) < p-f^'\ (29) 

for all r that satisfy 

r>„, I V (4^)'^' [bk + MT{1 + log p)f/' 
where is the lower k x k part of the upper triangular matrix R^^'^ of the QR decomposition of BiG, and /(r) = 



(M — r){N — r) for r £ [0, min{M, A^}). It must be noted that (30l is only valid for constant (fixed) bias term b, and a 



lower bound on F for the general expression of b (see ( 22 1 is known yet. However, as will be discuss in the sequel, even 
for fixed bias, the above lower bound shows a significant improvement in complexity compared to the more complex 
optimal sphere decoder. 

The above important result indicates that there exists a finite probability that the number of computations performed 
by the decoder may become excessive even at high SNR, irrespective to the channel being ill or well-conditioned! This 
probability is usually referred to as the probability of a decoding failure. Such probability limits the performance of the 
lattice sequential decoder, especially for a one-way communication system. For a two-way communication system, such 
as in our MIMO ARQ system, the feedback channel can be used to eliminate the decoding failure probability. Therefore, 
our new decoder must be carefully designed to predict in advance the occurrence of decoding failure to avoid wasting 
the time trying to decode the message. This would result in a huge saving in decoding complexity. As will be shown in 



the sequel, the above result can be easily extended to the MIMO ARQ channel. 

IV. TiME-OuT Algorithm 

It is well-known ||9l that the number of computations required to decode a message using sequential decoders is 
highly variable and assume very large values during intervals of high channel noise. As such, the decoder is expected 
to spend longer time attempting to decode the message. For the proposed incremental-redundancy MIMO ARQ system, 
this condition can be used as an indicator of when the receiver should terminate the search and request the transmitter 
for additional redundancy bits during any of the £ < L rounds. 

In order to avoid wasting time trying to decode a noisy signal during any of I < L ARQ rounds, we implement a 
time-out algorithm in the lattice stack sequential decoder for joint error detection and correction. Such algorithm works 
as follows: we define a parameter Tout to be the maximum time (number of computations) allowed to decode a message 
during any of the £ < L ARQ rounds. If the decoding time exceeds Tout, a NACK bit is fed back to the transmitter. 
The only exception of this rule is when the maximum number of ARQ rounds, L, is reached. In this case, the regular 
lattice sequential decoder (with no time-out limit) is used, where a NACK bit will be interpreted as an error, and the 
transmission of the next message is started anyway. Next, we define the retransmission probability and the undetected 
error probability from a lattice point of view. Those two quantities are responsible for the performance-throughput 
tradeoff achieved by the MIMO ARQ system. Throughout the work we assume the use of a small (fixed) bias term 
during all L ARQ rounds. Before continuing our analysis, we would like to introduce some important definitions related 
to the MIMO ARQ channel that will be used throughout the paper. 

Denote < A^''^ < • • • < ^li^in{M n] eigenvalues of {Hj)^Hj, VI < j < £. Let us first define the outage event 
for both long-term and short-term static channels under the lattice sequential decoder with bias term b that is given in 



(22 1 and r]{X, p) as defined in (20 1, with £ received blocks as 



Ob{p, £) = {m^ C^^*^ VI < i < ^ : Rb{p) < Ri 



where 

M 



Rh{p) 



£"^2 Ci log(l + /oAj), for long-term static channel; 



T\i (31) 
^ X] Ci(l + P-^i )) foTi' short-term static channel^ 



Denote, R\ = rilogp, and define Ui = — logAi/logp, and {x)~^ = max{0,x}., then It can be easily verified that at 
high SNR, the outage event for the long-term static channel model can be expressed as (assuming N > M) 



M ^ 

5]0(l-a.)+<ri/^l 
i=i ) 



(i) 

where in such channel we have A^- = Aj, VI < j < ^. For the short-term static channel we have 



I M 



• • • G Ri" : a? > • • • > ag, VI < j < A ^ E 



Ossii) = { 



Then, the associated asymptotic outage probabiUty is given by (see u5\i ) 



1 — a 



<ri} . 



j=i i=i 



Pout{p,i) 



Pi{Ois{i)) = p for long-term static channel; 

Pi{Oss{^)) = for short-term static channel, 
where dout(^) = db{ri/£) where db{r) is as defined in (?). 
A. Retransmission Request Probability 



(32) 



We make use of the lower and the upper bounds of the lattice sequential decoder's error performance in (26 1 to 



implement an error control mechanism for our IR-LAST MIMO ARQ system. It is clear from (13 1 and (26 1 that 
Vo{B£G) C Vo{B£G) for all b > 0. Therefore, at the decoder side, one can divide each of the Voronoi regions of the 
channel-code lattice A(B^G), i.e., Vu{B£G) (corresponds to a lattice point u = B^Gz, z G Z™) into two disjoint regions 
— TZuiB^G) and Vu{BeG)\Tlu{BeG). This is depicted in Fig. 1. For convenience, we define Vu{i) = VuiBeG) and 
Tluii) = TZuiBeG). 

Now, denote the region 'D{£) as 

V{i)=W^\} \J 7^nW>. (33) 

[ueA(BfG) J 

We take advantage of the feedback channel by introducing the erasure option at the decoder such that whenever the 
received signal yi € ^{i), for £ < L, the decoder requests for a repeat transmission or additional bits (as in the case of 
IR-LAST coding scheme). In this case, at round £ < L, the probability of a retransmission when the channel is not in 
outage is given by 

Pj:(Ai, no outage) = Pr(y£ G 'D{£)). (34) 

Unfortunately, determining whether G 'D{£) or not is considered by itself a difficult problem. We try to simplify this 
problem through the use of lattice sequential decoding by tracking the number of computations performed during the 
search for the closest lattice point. 



Following the definition of the number of computations performed by the decoder provided in (28 1, a retransmission is 
requested by the time-out algorithm at round £ < L if the number of computations exceeds the maximum time allowed 
before reaching the end of the tree. In other words, a NACK is sent back to the transmitter if, at any I < j < m, 
Mj{£) > Tout- This event could occur during a high channel noise period so that the received signal is close to the 
boundaries of a Voronoi cell of the lattice A{B£G) and as a result, the decoder declares that yi G ^(^)- In this case. 



for the selected value of b, one has to carefully choose the time-out parameter Tout so that whenever Mj{l) > Tont, the 
decoder decides that G ^(^) (see Fig. 2. (a)). Selecting an inappropriate value of Tout may result in the loss of the 
optimal tradeoff. Later, we shall make use of the following result for evaluating (34i, for fixecj^bias values: 



Lemma 1. For the long-term static ARQ channel, the asymptotic tail distribution of the total computational complexity 
of the lattice sequential decoder with fixed bias b > 0, at round i given the channel is not in outage, Pr(A/'j(£) > Tout). 
can be upper bounded by 

Pr(AA,(^) >rout) <P~*, Vl<j<m, (35) 



under the condition 



where C is a constant chosen sufficiently large enough so that MTLC, > (M — ri/i)(N — ri/£), 

AI 



inf i V(2i-1 + 7V-M)c 

M 



.i=l 



(36) 



and 



M 



Olsil) 



a G 



+ , ai > • • • > OA/ : ^(1 - ai)'^ > 



Proof: see Appendix L ■ 
It must be noted that the above result can be easily extended for the short-term static channel. In this case, one can 
show that Pr(7Vj(^) > Tout) < VI < j < m. Next, we derive an upper bound for the undetected error probability. 

B. Undetected Error Probability 

Another important parameter that characterizes the performance of the MIMO ARQ system, is the undetected error 
probability which was defined in ([8]). In the IR-LAST MIMO ARQ system, an ACK bit is sent back to the transmitter 
if the decoder correctly decode the received signal subject to the condition that Mm{(-) < Tout- This corresponds to the 
event G Tlo{i), assuming was transmitted (see Fig. 2.(b)). Also, an ACK is sent back to the transmitter, if decoding 
fails but it is not detected. Such an event occurs when the total number of computations Mm{(-) < Tout, but the decoded 
lattice point is not 0. This happens when the received signal yi G TZu{() for any u 7^ (see Fig. 2.(c)). Using the above 



The general case of a variable bias term that is given in \22\ will not be considered here due to the difficulty of obtaining a lower bound on 
Tout for such bias values. 



argument, the undetected error probability can be expressed as (assuming was transmitted) 

\ueA' {B,_G) J 



Our goal now is to upper bound ( |37| ). In this work, we will resort to a geometrical approach to obtain an upper bound 
on the undetected error probability. Before doing so, we express the event of sending an ACK, i.e., Ae, in terms of 
the number of computations performed by the lattice sequential decoder as Ae = {J\fm{i) < Pout}- Consider now the 
following theorem: 

Theorem 1. For any lattice code Ac, the undetected error probability of the quasi-static M x N MIMO ARQ system 
with maximum rounds L, and codeword length T, under the time-out MMSE-DFE lattice sequential decoding scheme 
with parameters h and Pout^ bounded above as 
In the case of long-term static channel 

Fi{Esd{Be,b),Ae)<p-f^'''/''\ (38) 
which is achieved with code block length T > \{M -\- N — 1)/L\, where f[r) = [M — r){N — r), and any b that 



satisfies {22). 



In the case of short-term static channel 

Pr(ii;,,(B,,6),^,)<p-^^(^'^/^), (39) 
which is achieved with code block length T > M -\- N — 1. 



Proof: It seems a bit difficult to obtain an upper bound directly for the undetected error probability using (37 1. 



Therefore, we resort to a geometrical approach (see Fig. 3) to further upper bound ( 37 1 by selecting b such that the 
bound VcsiBiG) < r^^c\{B(^G) is maintained for all £ < L. The effective radius VcsiBiG) of the lattice generated by 
BiG is given by 



Kdet(Bjs,)V2 



In this case, we can upper bound the undetected error probability as 



l/m 



(40) 



FiiEsdiBi,b),Ae) = Pr[ [j {e^ e 7l„(^)} ) < Pr(|e^|2 > r2g(S,G)). (41) 

^ueA* (BG) 



It is not yet clear how the RHS of (41 1 can be evaluated. To overcome this problem, we can find a lower bound on 
r'^g{BiG) at high SNR as follows: 



Asymptotically, one can express bias term that is defined by ( p2| ) as 



e/ML 



1 



V 



5Er=i(i-".)^ 



e/2ML' 



(42) 



Substituting (42) in (13 1, when the channelis not in outage, one can upper bound det(BjB^)) > r/^^, where 77 = p^^=l^'(^ 
In this case, we have that 



(a) 

> MTL 



Vc v(n) T 



2/m 



2/m 



> MTLp" > MTL{1 + 7 log p) , 



(43) 



where (a) follows from the fact that (see [2:3 1) there exists a shifted lattice code Ac + Uq with number of codewords 
inside the shaping region. 



|C(A„«S,7^)| = 2^^^ = p^^^> 



Also, 

^ — JTl l^^j—i ''^(''^ ~ '-^i)^ ~ ^1/^] ^ when the channel is not in outage, and the last inequality follows from 



the fact that limp^.oo(l + j log p)/p'^ = for any i^, 7 > 0. Therefore, 



Fr iEUBi,b),Ae) < Pr(|e^p > MLr(l + 7 log p)) < p-^^TL^. 



By choosing a large enough value of 7 such that MTL'y > /(ri/L), we obtain 



(45) 



(46) 



under the condition that LT > M + N - 1. 

The above analysis also applies to the short-term static channel, and one can show that under the condition T > 
N + M - I, we have 

PT{EUBi,b),A,) < p-^f(^'^/^). (47) 



As will be shown in the sequel, ( 34 1 and ( 37 1 play an important role in determining the achievable diversity- 



multiplexing-delay tradeoff of the MIMO ARQ channel. Reducing the number of retransmissions comes at the expense 
of increasing the undetected error probability, which is undesirable. It is clear that both probabilities are closely related 
and hence increasing or decreasing one of them may lead to a loss in the optimal tradeoff of the channel. 



C. Achieving the Optimal Tradeoff: Bias Term vs. Tout 

Our goal here is to prove the optimaUty of the time-out lattice sequential decoder in terms of the achievable diversity- 
multiplexing-delay tradeoff. It is well-known ll20l that the bias term b controls the amount of the computations performed 
by the lattice sequential decoder during the search and responsible for the excellent performance-complexity tradeoff 
achieved by such a decoder. Choosing a very large value of b although greatly reduces decoding complexity, it may 
lead to a loss in the optimal tradeoff of the channel. Therefore, it is expected that the time-out parameter Tout will be 
a function of the bias term b chosen in the algorithm. One may have already noticed that in Lemma 1. It turns out that 
an optimal value of Tout, denoted by r*^^^, exists so that the optimal tradeoff is achieved with a fairly low decoding 
complexity (i.e., average number of computations) compared to the joint list lattice decoder. The achievability of the 
optimal tradeoff, for any fixed bias term, under time-out lattice sequential decoding is summarized in the following 
theorem: 

Theorem 2. Consider a MIMO ARQ channel under short-power constraint given in (|2]), with M transmit, N receive 
antennas, a maximum number of ARQ rounds L, an effective multiplexing gain < rg < min{M, A^}. Then, the IR- 
LAST coding scheme under time-out lattice stack sequential decoding with parameter Tout <^nd fixed b > 0, achieves the 
optimal tradeoff: 



In the case of long-term static channel 



d*,{re,L) 



/(^), 0<re <min{M,iV}; 
0, re > min{M,7V}, 



which is achieved with code block length T > \(M -\- N — 1)/L\, where f{r) = (M — r){N — r). 
In the case of short-term static channel 



d*ss{''"e, L) 



( 

Lf[j), 0<re <min{M,iV}; 
0, re > min{M,7V}, 



which is achieved with code block length T > M + N — 1. The optimal tradeoff is achieved subject to the condition 

Proof: see Appendix III. ■ 
It must be noted that the above theorem is only valid for non-zero, but fixed values of b. Although one can fully 
characterize the achievable tradeoff in terms of the bias terrr|^(see Section III), for non-fixed bias, it is not yet clear how 



''For example, if we let 6 to vary with SNR and channel statistics as given in i22i one can achieve an asymptotic SNR exponent d;*^(re,L) 



(re/L) and desire, L) — Ldbife/L) for the long-term and short term static channels, respectively, where dh{r) is as described in section III. 



the time-out parameter Tout changes with the bias term for variable b. Therefore, in what follows, and for the purpose 
of completing the analysis, we will only consider the case of fixed bias term. 

Now, since Tout depends on the channel statistics (i.e., it is random), it would be desirable to determine its average 
value (averaged over channel statistics when it is not in outage). This may shed the light on determining the optimal 
value of Tout that can be used to achieve the optimal tradeoff of the channel. This can be done as follows: 

Consider the long-term static channel. Therefore, the optimal average value of Pout) say r*out> 

may be asymptotically 

lower bounded by 

m 

r^out = m + j;(logp)'=/2E,^o,.W {det(<f , (49) 

k=l 

In this paper we focus our analysis on nested LAST codes, specifically LAST codes that are generated using 
construction A that is described below (see |23l). 

We consider the Loeliger ensemble of mod-p lattices, where p is a prime. First, we generate the set of all lattices 
given by 

where p— S'oo,^— )-Oisa scaling coefficient chosen such that the fundamental volume Vf = ^^^^^p^^^^"^ = 1, 
Zp denotes the field of mod-p integers, and C C Z'^^'^'^^ is a linear code over Zp with generator matrix in systematic 
form [/ P'^]^. We use a pair of self-similar lattices for nesting. We take the shaping lattice to be = ^Ap, where ( is 
chosen such that the covering radius is 1/2 in order to satisfy the input power constraint. Finally, the coding lattice is 
obtained as Ac = />~''/^*^^As. Interestingly, one can construct a generator matrix of Ap as (see ll26l ) 

Gp = K , (50) 

\P plj 

which has a lower triangular form. In this case, one can express the generator matrix of Ac as G = where 
G' = (Gp. Thanks to the lower triangular format of G. If M is an m x m arbitrary full-rank matrix, where m = 2MTL, 
and G is an m X m lower triangular matrix, then one can easily show that 

det[(MG)fcfc] = det(Mfcfc) det(Gfcfc), (51) 

where {MG)kk, M^k, and Gkk, are the lower k x k part of MG, M, and G, respectively. 
Using the above result, one can express the determinant that appears in ([49]) as 



detiRi'^Ri'l) = det(B(f b2) det(Gj,Gfcfc) = p-^^ML a^^^R^jRW) det{G'l,G' ^k) , (52) 



(i)T (i) 

Let /^i < /i2 < • • • < ^fc be the ordered nonzero eigenvalues of Bj^j^ ^kk' for = 1, • • • , m. Then, 



det(<)Xi) = n/^. 

Note that for the special case when k = m we have fJ-2{j-i)TL+i = ••• = fJ'2jTL = 1 + p\j{{H^)^H'^), for all 
j = I,--- ,M. 



Denote a'- = — log /ij/ log p. Using (51 1, one can asymptotically express r*out as 



r*out = m + Y.i^ogpf'^E^^ou^M"}^ 



where 



Cfc 



k=l 



1 r + 

i=i 



(53) 



(54) 



Now, since is non-decreasing in k, we have 



(55) 



where 



M 
i=l 



At multiplexing gain ri, we have the channel is in outage only when Yl!jLi{^ — otj)^ < fi/L. In this case, we have 



{p'-} = j P'-Uoc) da 



-- {log pr^^ 
--{\ogpri''p'^'^\ 



exp logp 



M 



M 



1=1 



ML 



Oii 



1=1 



da 



where Osi{L) = {a G Rf : EfiiC^ " < n/^}, 



and 



Uri) = max 



M 



ML 



i=l 



It is not so difficult to see that the optimal channel coefficients that maximize (56 1 are 



(56) 



a* = 1, for i = 1, ■ ■ ■ , M — k, 



and 



Q*=0, for i = M- A; + !,••• ,M. 



Substituting a* in (56 1, we get 



for ri = 0, ,M. An since re = ri, the asymptotic average computational complexity, when the channel is in 
outage, can be expressed a^ 

T^out = 2MTL + (log . (58) 

One interesting special case of computing the optimal average time-out parameter r*out is when rg = 0, i.e., when 
using a code with fixed rate Ri. In this case we have 

r*out = 2MTL + ^ . (59) 

Thus, the above equation describes how the time-out limit is related to the system parameters (M, N, T, L, p) at the high 
SNR regime. As an example, consider a 2 x 2 MIMO ARQ channel with L = 2 rounds. Then, according to Theorem 
2, we have T = 3 is sufficient to achieve the optimal tradeoff. The signal dimension in this case is m = 24. Assume 



p = 100 (20 dB). According to (58 1, the optimal time-out limit is given by r*out > 98. As will be shown in the 
sequel, this theoretical result closely matches the value of r*out that is obtained experimentally (see Section V). Typical 
values of b that corresponds to r*out ~ 98 are between 0.6 and 1. For very small values of b, the average number of 
computations increases and according to that r*out ^ 98. 

To see the great advantage of using the time-out lattice sequential decoder with constant bias term over the list 
lattice decoder implemented via sphere decoding algorithms, we compare the average computational complexity of both 
decoders. It has been shown in 1127 II that, for moderate-to-high SNR, the average computations performed by the sphere 
decoder when the channel is not in outage, say Fsphere for a system with m = 2MTL signal dimension is given by 
(assuming fixed rate ri = 0) 

Tsphere = 2MTL + ^^J^ • (60) 

The ratio of the average complexity of both decoder, say 7, is given by 



r^phere ^ 2MTL + (log p)^^^^^/^^ 

r^„,t 2MTL + (log pY^TL/pMN ■ 



This is a huge saving in computational complexity, especially for large signal dimensions even at very high SNR. For 



^As a reminder, the logarithm that appears in the complexity analysis is to the base 2. 



example, at p = 10 (80 dB), 7 « 7.4, i.e., the list sphere decoder's complexity is about 7 times larger than the 
complexity of the new proposed time-out lattice sequential decoder. For p < 80 dB, one would expect the ratio 7^7. 
For extremely high SNR values (e.g., p > 90 dB), it seems that 7 — 1 as p — 00. 

Moreover, some interesting remarks about the effect of extremely decreasing or increasing the value of b on the 
performance-throughput-complexity tradeoff are discussed next. For very large SNR values, it is 

It is well-known that as fe — )• 00, lattice sequential decoders based on Schnnor-Euchner enumeration converts to 
the MMSE-DFE decoder EOl . In terms of the total number of visited nodes, MMSE-DFE decoder achieves linear 
computational complexity in m. In this case, for any Tout > "i, the message will always be decoded from the first 
round. Although it achieves high throughput and is computationally efficient, this decoder cannot achieve the optimal 
tradeoff. Assuming N > M, the maximum SNR exponent that such a decoder can achieve is {N — M + l){l — re/ML)^ 
(see [22 1 for more details about MMSE decoding for the case of MIMO channel with no ARQ). 

On the other hand, as 6 — )• the decoder achieves the best performance. However, the decoding complexity becomes 
equivalent to (and for some cases worst than) the complexity of lattice (sphere) decoding which is extensively large. Our 
main objective of using lattice sequential decoding is to save on decoding complexity. Therefore, one should appropriately 
select the bias term b so that the optimal tradeoff is achieved while reducing the computational complexity. This may be 
achieved by the sequential decoder by ensuring that the path metric along the correct path increases on average, while 
decreases along other paths. In this case, we choose b such that ¥,gr{fj.{zi)} > (assuming is the correct path). This 
corresponds to 6 > E{|[e']jp} = 1/2. This fact is verified experimentally as will be shown in the next section. 

As discussed earlier, the value of the parameter Tout used in the time-out algorithm is critical for achieving the optimal 
tradeoff of the MIMO ARQ system. In order to achieve the optimal tradeoff, both b and Tout have to be appropriately 
selected in the time-out algorithm. Remember that the probability of error is upper bounded by 

L-l 

Pe<Yl y^iEe,Ae) + PrCEO (61) 
controlled by ro„t controlled by b 

where PT{E£,Ae) < p~-f(^'/^') is the probability of undetected error at round £ and is mainly controlled by the time-out 
parameter Tout- However, the second term of the RHS of the above upper bound is affected by the value of the bias term 
b. Therefore, both b and Tout have to be appropriately selected in the time-out algorithm so that a balance is obtained 
which may lead to achieving the optimal tradeoff of the channel. 

Suppose that an optimal value of b is set in the time-out algorithm. Then, choosing a very small or large value of Tout 
may result in a loss of the optimal tradeoff. This is because Tout can be seen as the parameter responsible for the amount 
of retransmission when error is detected. Choosing a large value of Tout reduces the probability of retransmission and 
may result in large performance degradation. In this case, as Tout — ^ 00, we have Pr(^^) — )• 1. Therefore, the undetected 



error probability becomes equivalent to (for the long-term static channel) 

Vr{Ei,A,)=VT{E,)=p-f^'-^l'\ 



Moreover, the probability of retransmission given the channel is not in outage approaches 0, which means that = ri 
{p{tj — )• 0). And since the average error probability defined in ([S]) is dominated by the term with the smallest SNR 
exponent, we have 

Pe{p) < + ^p-fir^m = p-fir^). (62) 

e=i 

This is equivalent to the performance of the MIMO channel with no ARQ. On the other hand, choosing a very small 
value of Tout improves the performanc^ at the expense of large throughput loss. In this case, the undetected error 
probability approaches 0, p{£) — )• 1, and the throughput rj — )• Ri/L, i.e., Ve = ri/L. The error probability at high SNR 
in this case is also given by ([62]). 

Therefore, in order to achieve the optimal tradeoff of the MIMO ARQ channel, first b must be chosen to ensure that 
the achievability of the optimal diversity-multiplexing tradeoff when operating over the whole received signal (i.e., at 
round L). Then, Tout is selected accordingly so that the optimal diversity-multiplexing-delay tradeoff is achieved. 

V. Simulation Results 

In our simulation, we consider a long-term static MIMO ARQ link with M = N = L = 2, T = 3 and i?i = 8 

bits per channel use. The incremental redundancy LAST code is obtained as an {m,p,k) Loeliger construction (refer 

to ll23l for a detailed description). The frame error rate and computational complexity are plotted in Fig. 4 and Fig. 5, 

respectively, for different values of b used in the time-out algorithm. We measure the computational complexity of the 

joint lattice sequential decoder as the total number of nodes visited by the decoder during the search, accumulated 

over all ARQ rounds, until a new transmission is started. For every value of b, the optimal value of Tout, denoted 

by T*^^, was found via simulation by trial and error Using those optimal parameters, it is shown that the IR-LAST 

coding scheme decoded using the time-out lattice sequential decoder can achieve probability of error very close to the 

one that corresponds to the same IR-LAST coding scheme decoded using the list lattice decoder This is achieved with 

significant reduction in complexity compared to the list lattice decoder (see Fig. 5). It is interesting to see how such 

IR-LAST coding scheme can achieve probability of error colose to the coherent LAST code with half the rate (4 bpcu). 

On the other hand, the effective rate, i?e, of the IR-LAST coding scheme decoded under the new proposed decoder is 

shown to approach i?i = 8 as SNR grows as predicted by the theory. Optimal values of Tout for some special cases of 

b are provided in Table 1. As expected, for values of 6 < 1/2, the average computational complexity of the time-out 

^Performance improvement is acliieved as a coding gain and not in the SNR exponent. In tliis case, the decoder will accumulate more information 
about the message before making a decision which is due to the fact that most of the time the decoder asks for retransmission. 



algorithm increases and as a consequence, the value of r*^^ is proportionally increased. Simulation results demonstrate 
the excellent performance-complexity tradeoff achieved by the proposed algorithm for all values of b, especially at the 
moderate-to-high SNR regime (see Fig. 5). 

TABLE I 

Optimum values of Tout for some special cases of b used in the time-out algorithm for the case opM = N = L = 2 and 

T = 3 MIMO ARQ system USING IR-LAST random code 



b 


^ out 


0.6 


100 


0.4 


800 


0.1 


4 X 10^ 



The error rates are obtained by averaging over at least 10 000 channel realizations at small SNRs and as much channel 
realization as required to count at least 100 frame errors at high SNRs. It is clear that for values of 6 > 1/2 the decoder 
exhibits some performance degradation compared to the incomplete list sphere decoder. To improve the performance 

one need to resort to smaller values of b at the price of increasing computational complexity (see Fig. 4 and Fig. 5). 
In general and for all finite values of b, the time-out lattice stack sequential decoder has much lower computational 
complexity compared to the incomplete list sphere decoder. 

VI. Conclusion 

In this paper, we have demonstrated, analytically and via simulation, the significant improvements achieved by the 
lattice stack sequential decoder over the incomplete hst lattice decoder that is used for joint error detection and correction 
in the IR-LAST MIMO ARQ channel. A time-out algorithm has been proposed. Theoretical analysis and simulation 
results show that the optimal tradeoff can be achieved by such algorithm with very low decoding complexity compared 
to the hst lattice decoder, especially for moderate-to-high SNR for which the list output could be extensively large. 

Appendix I: 
Proof of Lemma 1 

Let <p'{zi,£) be the indicator function defined by 

1, if |e^'{-i224l'<^^-/^min(^); 

0, otherwise, 

where /Uinin(^) is the minimum metric along the decoded path. Then, it can be easily verified that 



where (t){zi,£) is as defined in (27i, where a path may be extended by the stack decoder if the partial path metric at 



that node satisfies > Mmiii(^)- Now, given \e'^\^ < MTL{1 + 7), and by noticing that -(/UminW + ) < 0, 



we obtain 



where 



2f C 



1, if \ef^ - Rtlz\? <bk + MTL{1 + 7); 
0, otherwise. 



Notice the independence of the upper bound (64i on /imin(^)- Now, let 



Sk{i), if \e'l - Rez\^ < bm - fi^U^y, 
0, otherwise. 



where 



then it can be easily shown that 



Sk{i)= Y '^"(^i' 



2*ez'= zeZ" xeA, 



where 



S'fc(£), if |S£a;p - 2(S^x)"^e^ < 6m; 
0, otherwise. 



Consider now the following lemma: 



(64) 



(65) 



(66) 



Lemma 2. At the i-th ARQ round, the number of nodes visited by the lattice stack sequential decoder at level k 
1 , • • • , m, given that |e'p < AfrL(l +7), can be upper bounded by (for any finite b > 0) 



(4„f''' [bk + MTHl + 7)]'=/2 



ik/2 + l) det(i2ifi?W)V^ 



(67) 



where Sk{i) is as defined in (661. 



Proof: See Appendix II. 
The tail distribution, given the channel is not in outage, can then be upper bounded as follows 



Pr(ACW > ToutPUl)) < Pr(|e^|2 > MTL{1 + 7)) + Pr(A/'„(^) > Tout, |e^P < MTL{1 + ^)\OUe)). (68) 



For a given lattice Ac, using Markov inequality, we have 



Pr(AA„(^) > ToutlAc, Oisit), |e^r < MTL{1 + 7)) < Pr(AA„(^) > Tout - m|Ae, |e^|' < MrL(l + 7)) 

Ee.{AC(^)|Ac,OisW,KP <MrL(l + 7)} 

< , for Tout > m. 

Tout - m 

(69) 

where Mm{^), assuming all- zero codeword was transmitted, is defined as 

m 

fc=iztez''\{0} 

Using Lemma 2, the conditional average of Mm{^) with respect to the noise can be further upper bounded as 

m 

Ee'{M„(£)|Ae, |e^p < MTL{1 + 7)} < J] J] Pr(|B^x|2 - 2{Bix)^e'^ < bm). (70) 

k=l xeA* 

Therefore, 

m 

Pr{Mmie) > ToutlAe, O/sW, |e^|' < MrL(l + 7)) < ^ V Fri\Bex\^ - 2(B£x)V^ < bm). (71) 

Tout - m ^ 

We would like now to upper bound the term inside the summation in ( [TT] ). The difficulty here stems from the non- 
Gaussianity of the random vector for any finite T. To overcome this problem, consider the following: 
Let 

ee = [Bi - FiHi]gi + Fe{ei + Wi), 

where ge ~ M{0,a'^lm), Wi ~ M{0, (a^ — 1/2)/^) and cj^ > 1/2. Following the footsteps of lHH, it can be shown 
that by appropriately constructing a nested LAST code we have that 

Pr(|B^a;p - 2(S<.x)"^e^ < bm) < Pr(|B£a;p - 2(S^x)"^e£ < bm), (72) 

where ~ A/'(0, 0.51m), and f3m is a constant independent of p. Using Chernoff bound, 

{g-|S.rr|V8gfem/4^ \B^x\^ > 6rn; 
(73) 
1, \Bgx\^ < bm. 

By taking the expectation over the ensemble of random lattices (see flS], Theorem 4) 

m/2(t^ \m/2 to \m/2 bm/i ^ 



Substituting ^ and ^ in and by selecting Tout >m + Yl'k=i S'ki^), we get 



(75) 



Now, by setting 7 = ^logp, the first term in the RHS of (681 can be shown to be upper bounded by p /(^i/^), as long 
as ( is chosen sufficiently large such that MTLQ > f{ri/£) (see |[T5l . Appendix IV). Averaging the second term in the 



RHS of ( 75 1 over the channels in Ois {£) set, we obtain , 



PriMmii) > Tout) < p-^^'''^'^ + [_ Ua) Pr(A/-„ > Toutl") <i« < p-f^''"\ (76) 
where fa{ct) is the joint probability density function of a which, for all a G Ois{l), is asymptotically given by (see fSM ) 

_ (77) 

i=l / 

Since Nj{() < Nm{t), for all I < j < m, then 



Ua) = exp ^- log{p) Y,{2i -l + N- M)a}j 



PiiMjil) > Tout) < Pr(AA„(^) > Tout) < p-^^''/'\ 

Appendix II: 
Proof of Lemma 2 

In the proof of Lemma 1, we have shown that computational complexity at level k when |e^p < MTL{1 + 7), can 
be upper bounded by 

where (f)"{zi,£) is the indicator function defined in (65 1. For brevity, we drop the subscript £ in the rest of the proof. 



Given that |ep < ii^, it must follow that {e^l < R^, where e'f is the last k components of e. Without loss of generality, 
we assume that all-zero lattice point was transmitted. Let 



(4) = 



1, if |e't - Rkkz'tl^ <bk + Rl |e't|2 < i?2; 
0, otherwise. 



(78) 



where R^ = MTL{1 + 7). The total number of integer lattice points that satisfy (78) is given by 



(79) 



where 

,1, if \e't - Rkkz'tl^ <bk + Rl\e't\^ <bk + R^,; 
'"{z'l) = { (80) 

0, otherwise. 

In general one can show that for any random vectors u and v, and a > 0, it holds{|tt— < a, \u\'^ < a} CI {\u\'^ < 4a}. 



Therefore, we can further upper bound ( [79] ) as 
where 



{z\) 



1, if\RkkzlV <A{hk + Ri); 

(82) 



0, otherwise. 

The summation of (p{zi) over all integer lattice points z\ G Z'^ can then be easily upper bounded by (see [13]) 



Appendix III: 
Proof of Theorem 2 

First, the error probability (|8]l is lower bounded by the probability of error of the optimal maximum-likelihood or the 
MMSE-DFE lattice decoder that operates on the whole received signal vector y = Ul knowing the channel matrix H'^. 
Hence, it can be easily shown that (see ifTOl or O) 

Pe{p) > P{Eid{BLG)) = Pout(p, L). (83) 

The input to the decoder at round I, after QR preprocessing of (|4]) is given by 

ye=Qjye=Rez + el (84) 

where e" = Qje'^. 

Consider the long-term static channel. For the probability of error we bound each term in ([8]l. At round L, due to 
lattice symmetry, we assume that the all zero codeword, i.e., 0, was transmitted. At high SNR, we have shown in section 
III that, for a given lattice Ac, the error probability of the lattice stack sequential decoder can be upper bounded as 

P{EUBl, b)\A,) < P{Eia{BL)\A,). (85) 



where Bl = {I — €)Bl, and < e < 1. Using the union bound, we can further upper bound (85 1 by 



P{Esd{BLM^c) < Pr(2(SLx)V > \Blx\ 



(86) 



As shown in the proof of Theorem 1 (see Appendix I), by appropriately constructing a nested LAST code we have 
that 

P{EsdiBL, 6)|Ae) </3mYl Pr(2(^Lx)Te > \Blx\^), (87) 



xeA* 



(88) 



where e ~ A/'(0, 1/21^), and /3m is a constant independent of p. Using Chernoff bound, 

Fii2{B'Lxfe > \Blx\^) < e-l-^-^l'/^. 

By taking the expectation over the ensemble average of random lattices (see {23], Theorem 4) 

P{E,d{BL, b)) = EAAP{Esd{BL, b)\Ac)} 

- V, J - Vcil-e)--/^detiBlBLy/^' 

\B'^x\^>0 

Substituting ^ (with £ = L) and (44i in (89 1, the average error probability at round L, given the channel is not in 
outage can be rewritten as 

P{E,d{BL,b)\OUL)) < /C(m,e)p-^^E;-(i-"^)"-xl, (90) 



(89) 



where /C(m,e) is a constant independent of p. Using the results in lITOl . the average probability of error (average over 
the channel and lattice ensemble) is asymptotically upper bounded by 



P{E,d{BL, b)) < FiiOUL)) + P{Esd{BL, b),Ou{L)) 



< 



P 



+ 



[_ Ua)P{Esd{BLMOis{L)) dn < p-f^'^/'^l 

JOUL) 



under the condition that LT > AI + N — 1, where /q.(q;) is given in (77 1. 

The proof for the short-term static channel follows the same lines with the exception that 

L 



det{BlBL) = det (/ + {^{H'^jfH^j 



2T 



i=i 



In this case, one can easily verify that 



(91) 



(92) 



under the condition that T > M + N - I. 



Thus, the error probabiUty in Q can be asymptotically upper bounded as 

/or long-term static channel] 
p-Lf{ri/L) ^ /or short-term static channel. 



Peip) < 



Next, we need to show that = ri. As discussed previously, when the channel is not in outage, a retransmission is 
requested by the time-out lattice sequential decoder at round £, whenever the number of computations Mj (i) at any level 
1 < j < m exceeds Tout- Therefore, the probability of a retransmission at round £ when the channel is not in outage, 
can be expressed as 

And since Nm{t) > J^j{£) VI < j < m, using the union bound we can upper bound the above probability by 

Y>i{At,Ois{pA)) < mY>^{Mm{£) > Tout). 

Therefore, we have 

p{£) < Ft(A,) 

= FT{OUp,£)) + FT(Ai,OUp,e)). (93) 
The behavior of the first term at high SNR is p~^^'^'^/^\ Using Lemma 1, we get 

VT{A,,Ois{p,£))<p-'''. 

where is as defined in Theorem 1. 

Following 1 10], for any T > 1 we find that > for all £ and n < M. Moreover, if T£ > M + N - I then 
di = f{'ri/£), which is the maximum possible SNR exponent for codes with multiplexing gain ri/£ and block length 
T£. 



Therefore, we have that 



Using ([T]), we obtain 



1 + Eti p 



'mm{f{ri/e),di} 



<ri<Ri. 



This implies that = ri, and therefore for the long-term static channel, (i^^(re,L) is achievable. The proofs for the 
short-term static channel follow exactly the same arguments, and one can show that d*g{re,L) is achievable under 



time-out lattice sequential decoding. 
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Fig. 1. The division of the Voronoi cell of the lattice generated by the channel-code matrix B(G into two distinct regions — the shaded region 
TZu{BeG), and the white region Vu{BeG)\TZu{BeG). The two dimensional hexagonal lattice is shown for illustration purposes. 




(a) The event of sending a NACK when 
the channel is not in outage. 




(b) The event of sending an ACK with (c) The event of sending an ACK with 
correct decoding. decoding failure but not detected. 



Fig. 2. The events of retransmission, correct decoding, and undetected error that occur in the time-out algorithm (assuming was transmitted). 
The correct decoding region is represented in dark color. The chessboard shaded regions represent the undetected error events {E£,Ai). The white 
region represents the undetected error event. 




Fig. 3. A geometric approach used to over bound the undetected error probabiUty under the time-out algorithm. 
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Fig. 4. The optimal tradeoff achieved by the time-out algorithm lattice stack sequential decoder for several values of h. 
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Fig. 5. Comparison of average computational complexity of the MMSE-DFE list lattice decoder and the time-out lattice stack sequential decoder 
for several values of b using their corresponding optimal values of Tout (see Fig. 4). 



